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ABSTRACT 


The  shape  of  a  two-dimensional  bubble  rising  at  a  constant  velocity  U 
in  a  tube  of  width  h  is  computed.  The  flow  is  assumed  to  be  inviscid  and 

incompressible.  The  problem  is  solved  numerically  by  collocation.  The 

,  ' . \ 

results  confirm  Garabedian's  [2]  findings.  There  exists  a  unique  solution  for 

1/ 

each  value  of  the  Froude  number  F  =  U/(gh)£f,  smaller  than  a  critical  value 

y  t  f  i  *  '>!*.*>  ■ — 

~E_^_  Here  g  denotes  the  acceleration  of  gravity.  It  is  found  that  = 
0.36.  In  addition  the  problem  of  a  jet  emerging  from  a  vertical  nozzle  is 


considered.  It  is  shown  that  the  slope  of  the  free  surface  at  the  separation 

c.  c. 

points  is  horizontal  for  F  <  F^  and  vertical  for  F  >  F^.  Graphs  and  tables 
of  the  results  are  included.  __ 
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SIGNIFICANCE  AND  EXPLANATION 


We  consider  the  flow  of  an  incompressible,  inviscid,  heavy  liquid  past  of 
a  gas  bubble  in  an  infinitely  long  vertical  tube  (see  Figure  la).  The  flow  is 
assumed  to  be  two-dimensional.  This  problem  was  first  considered  by  Birkhoff 
and  Carter ^ ^  and  Garabedian^.  Birkhoff  and  Carter^ 1 ^  suggested  that  the 
problem  had  a  unique  solution.  On  the  other  hand  Garabedian'  '  presented 
analytical  evidence  that  the  problem  had  a  solution  for  each  value  of  the 
Froude  number  F  »  U/(gh)^2  smaller  than  a  critical  value  Fc.  Here  U  is 
the  velocity  at  infinity,  g  the  acceleration  of  gravity  and  h  the  width  of 
the  tube. 

In  the  present  paper  we  settle  the  controversy  between  the  results  of 

Birkhoff  and  Carter  ^ 1 ^  and  those  of  Garabedian  ^ 2  * .  We  compute  accurate 

(2) 

solutions  by  a  collocation  method.  Our  results  confirm  Garabedian 1 s ' 
findings.  There  exists  a  unique  solution  for  each  value  of  F  smaller  than  a 
critical  value  Fc.  However  we  found  Fc  *  0.36.  This  value  is  about  40 
percent  higher  than  that  indicated  by  earlier  work  on  the  problem. 

In  addition  we  consider  the  problem  of  a  jet  emerging  from  a  vertical 
nosxle  (see  Figures  1b  and  1c).  We  show  that  the  slope  of  the  free  surface  at 
the  separation  points  is  horizontal  for  F  <  Fc  and  vertical  for  F  >  Fc. 


The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive 
summary  lies  with  MRC,  and  not  with  the  author  of  this  report. 


BUBBLES  RISING  IN  A  TUBE  AND  JETS  FALLING  FROM  A  NOZZLE 


Jean-Marc  Vanden-Broeck 


I«  Introduction 

We  consider  the  two-dimensional  flow  of  an  incompressible  fluid  past  a 
gas  bubble  in  an  infinitely  long  tube  of  width  h.  We  assume  that  the  bubble 
extends  downwards  without  limit.  We  choose  a  frame  of  reference  moving  with 
the  bubble,  so  that  the  flow  at  large  distances  from  the  bubble  is 
characterized  by  a  constant  velocity  U  [see  Fig.  (la)].  We  choose  the 
origin  of  coordinates  at  the  top  of  the  bubble.  We  assume  that  the  flow  is 
symmetric  about  the  x-axis  and  that  gravity  g  is  acting  in  the  negative  x- 
direction.  As  we  shall  see  that  the  shape  of  the  bubble  is  determined  by  the 
Froude  number 


F 


(1) 


1  2 

This  problem  was  first  considered  by  Birkhoff  and  Carter  and  by  Garabedian  . 

Birkhoff  and  Carter1  attempted  to  solve  the  problem  numerically.  They 
assumed  that  the  problem  had  a  unique  solution  and  that  the  Froude  number  F 
could  be  obtained  as  part  of  the  solution.  Though  they  obtained  approximate 
solutions  with  F  ~  0.23,  the  convergence  of  their  procedure  was  not  really 
satisfactory. 

2 

Garabedian  presented  analytical  evidence  that  the  solution  is  not 
unique.  He  suggested  that  a  solution  exists  for  each  value  of  F  smaller 
than  a  critical  value  Fc.  In  addition  he  showed  that  Fc  >  0.2363  and 
guessed  the  value  Fc  *  0.24. 


Sponsored  by  the  United  States  Army  under  Contract  No.  DAAG29-80-C-004 1 . 
Also,  by  the  Australian  Research  Grant  Committee  and  the  National  Science 
Foundation  under  Grant  No.  MCS-8001960. 


In  the  present  paper  we  confute  accurate  numerical  solutions  by 

collocation.  Our  scheme  is  an  improved  version  of  the  procedures  proposed  by 

Birkhoff  and  Carter  <  Our  results  confirm  Garabedian's  findings.  There 

exists  a  unique  solution  for  each  value  F  smaller  than  a  critical  value 

F_.  However  we  found  that  F_  =  0.36.  This  value  is  about  40  percent  higher 
c  c 

2 

than  the  value  guessed  by  Garabedian  . 

The  flow  configuration  of  Fig.  1(a)  can  also  serve  to  model  a  jet 


rging  from  a  nozzle  [see  Fig.  1(b)].  As  x  tends  to  infinity,  the 


velocity  is  assumed  to  approach  a  constant  U.  It  follows  from  the  symmetry 
of  the  flow  configurations  that  the  portion  SJ  of  the  bubble  surface  in 
Figure  la  is  identical,  for  the  same  value  of  F,  to  the  portion  SJ  of  the 
jet  surface  in  Figure  1b.  Therefore  our  results  for  the  rising  bubble  imply 
that  the  jet  of  Figure  1b  exists  for  all  values  of  F  <  Fc.  For  F  >  Fc, 
this  flow  configuration  fails  to  exist. 

For  F  >  Fc  we  sought  a  solution  in  which  the  flow  separates 
tangentically  from  the  nozzle  [see  Fig.  1(c)].  We  found  that  there  exists  a 

unique  solution  of  the  type  sketched  in  Fig.  1(c)]  for  each  value  of 

F  >  F_.  For  F  <  F_  these  solutions  fail  to  exist. 

c  c 

From  our  results  we  conclude  that  there  exists  a  unique  jet  for  each 
value  of  F.  For  F  >  Fc  the  slope  of  the  free  surface  is  vertical  at  the 

separation  points  S  and  S'  For  F  <  Fc  the  slope  of  the  free  surface  is 

horizontal  at  the  separation  points  and  the  velocity  at  these  points  is  equal 
to  zero. 

The  problem  is  formulated  in  Sec.  II.  The  numerical  procedure  is 


described  in  Sec.  Ill  and  the  results  are  discussed  in  Sec.  IV. 
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IZ.  Formulation 


Let  us  consider  a  two-dimensional  jet  emerging  from  a  nozzle  of  width 
h  [see  Figs.  1(b)  and  1(c)].  As  x  •,  the  velocity  approaches  a 
constant  U.  The  fluid  is  assumed  to  be  incompressible  and  inviscid. 

We  define  dimensionless  variables  by  choosing  h  as  the  unit  length 
and  U  as  the  unit  velocity.  We  introduce  the  potential  function  t  and  the 
stream  function  4>.  Without  loss  of  generality  we  choose  $  *  0  at 
x  *  y  ■  0  and  $  =  0  on  the  streamline  IJ.  It  follows  from  the  choice  of 
the  dimensionless  variables  that  i|>  =  -1  on  IJ'.  The  complex  potential 
plane  is  sketched  in  Fig.  2. 

We  denote  the  complex  velocity  by  5  =  u  -  iv  and  we  define  the  function 
x  -  i0  by 

C  =  u  -  iv  =  eT  .  (2) 

We  shall  seek  T  -  i0  as  an  analytic  function  of  f  =  $  +  i+  in  the  strip 
-1  <  <  0. 


On  the  surface  of  the  jet,  the  Bernoulli  equation  yields 


1  2 


1  2 


2  q  +  **  "  2  qc 


(3) 


Here  q  is  the  flow  speed,  g  the  acceleration  of  gravity  and  qc  the 

velocity  at  the  separation  points.  In  dimensionless  variables  (3)  becomes 

2 

(4) 


2t  2 

e  +  —  x  =  —  on  SJ  and  S'J' 
FZ  U2 


Here  F  is  the  Froude  number  defined  by  (1). 

It  is  convenient  to  eliminate  y  and  qc  from  (4)  by  differentiating 
(4)  with  respect  to  Using  the  relation 


+  i  &  =  _ 1 
A  1  A  A  1  _ 


-T+i0 


3$  3$  u-iv 


(5) 


we  obtain 


e2T  tt  +  e  T  cos  8  “  0  on  SJ  and  S'J'  . 
*  F 


(6) 
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This  function  must  be  analytic  in  the  strip  -1  <  i>  <  0  and  satisfy  the 
conditions  (6)  -  (8)  and  (9)  for  the  flow  configuration  of  Fig.  1(b)  and  the 
conditions  (6)  -  (8)  and  (10)  for  the  flow  configuration  of  Fig.  1(c). 


III.  Numerical  procedure 

Following  Birkhoff  and  Carter1  we  define  the  new  variable  t  by  the 
relation 


-nf  _  1  ...  1. 
e  ■  y  (t  +  — ) 


(11) 


This  transformation  maps  the  flow  domain  onto  the  unit  circle  in  the  con$>lex 
t-plane  so  that  the  walls  of  the  nozzle  go  onto  the  real  diameter  and  the  free 
surfaces  onto  the  circumference  (see  Fig.  3). 

In  order  to  obtain  solutions  for  the  flow  configuration  of  Fig.  1(b),  we 
note  that 

1/3 


C  ~  [An(1  ±  it)}  '  as  t  ±i  (12) 

C~1±t  as  t+^1  (13) 

(see  Birkhoff  and  Carter1  for  details).  Therefore  we  define  the  function 
fl(t)  by  the  relation 

eT"i6  -  -(-An  C( 1  +  t2)]1/3(-An  C)‘1/3  (1  -  t2)en(t)  .  C 14) 
Here  C  is  an  arbitrary  constant  between  0  and  0.5.  We  choose  C  *  0.2. 
The  function  fl(t)  is  bounded  and  continuous  on  the  unit  circle  and  analytic 
in  the  interior.  The  conditions  (7)  and  (8)  show  iJiat  ft(t)  can  be  expanded 
in  the  form  of  a  Taylor  expansion  in  even  powers  of  t.  Hence 


00 

e"1"1  -  -[-An  C(1  +  t2)]V3(-An  C)”1/3  (1-t2)exp(£  a  t2n)  .  (15) 


n=1 


The  functions  t  and  0  defined  by  (15)  satisfy  (7)  -  (9).  The 
coefficients  an  have  to  be  determined  to  satisfy  (6)  on  SJ.  The  condition 
(6)  on  S’J’  will  then  be  autotomically  satisfied  by  symmetry. 


We  use  the  notation  t  ■  I t | e^a  so  that  points  on  SJ  are  given  by 

IT 

T 


t  **  e*0 ,  0  <  o  <  ■£•.  Using  (11)  we  rewrite  (6)  in  the  form 


_  .  „  2t  dr  1  -t  ~ 

x  cotg  0  e  - - —  e  cos  0  =  0 


(16) 


Here  t(o)  and  0(0)  denote  the  values  of  t  and  9  on  the  free  surface 

SJ. 

We  solve  the  problem  approximately  by  truncating  the  infinite  series  in 
(15)  after  N  terms.  We  find  the  N  coefficients  an  by  collocation.  Thus 
we  introduce  the  N  mesh  points 

°I  "  (I  ”  5)f  1  “  1  .  (17) 

iw 

Using  (15)  and  (17)  we  obtain  [t(o)]  ,  [0(c)]  and  I~)  in  terms 

0=0^  o=  Cj  do  o*o ^ 

of  the  coefficients  an<  Substituting  these  expressions  into  (16)  we  obtain 
N  nonlinear  equations  for  the  N  unknowns  an,  n  *  1,...,N.  We  solve  this 
system  by  Newton's  method.  Once  this  system  is  solved  for  a  given  value  of 
F,  we  calculate  the  functions  T  and  9  by  setting  |t|  *  1  in  (14).  The 
shape  of  the  jet  is  then  found  by  numerically  integrating  the  relation  (5). 

We  shall  refer  to  this  numerical  procedure  as  scheme  I. 

Solutions  for  the  flow  configuration  of  Fig.  1(c)  can  be  obtained  by 
omitting  the  factor  1  -  t  in  (14).  Thus  (15)  becomes 

eT_10  -  -[-fn  C(  1  +  t2)JV3(-*n  C)_1/3  exp(  \  b  t2")  .  (18) 

n=1  n 

The  functions  T  and  9  defined  by  (18)  satisfy  (7),  (8)  and  (9).  The 
coefficients  bn  have  to  be  found  to  satisfy  (6)  on  SJ.  We  truncate  the 
infinite  series  in  (18)  after  N  terms  and  determine  bn,  n  *  1,...,N  by  the 
collocation  method  described  in  scheme  Z.  We  shall  refer  to  this  numerical 
procedure  as  scheme  II . 

Solutions  for  the  flow  configuration  of  Figs.  1(b)  and  1(c)  can  also  be 
obtained  by  assuming  the  expansion 

eT’i0  -  -l-*n  C(  1  +  t2)]  1/3  (-An  cfV3  (1  +  \  dn  t2n)  . 

n=1 


(19) 


It  can  easily  be  verified  that  (19)  satisfies  (7)  and  (8).  In  addition  (19) 

«  m 

satisfies  (9)  when  7  d  -  -1  and  (10)  when  7  d  +  -1.  Therefore  (19) 

n“1  n*0 

is  appropriate  to  describe  the  flow  configuration  of  Fig.  Kb)  as  well  as  the 
flow  configuration  of  Fig.  1(c).  The  infinite  series  in  (19)  is  truncated 
after  N  terms  and  the  coefficients  d^,  n  ■  1,...,N  are  determined  by  the 
collocation  method  described  in  scheme  I.  We  shall  refer  to  this  scheme  as 
scheme  III. 


We  used  scheme  I  to  compute  solutions  for  the  flow  configuration  of  Fig. 
1(b).  We  obtained  accurate  solutions  for  F  <  0.3.  However  we  were  unable  to 
obtain  reliable  solutions  for  F  >  0.3.  Similarly  we  obtained  accurate 
solutions  for  the  flow  configuration  of  Fig.  1(c)  for  F  >  0.4  by  using 
scheme  II.  However  we  could  not  calculate  reliable  solutions  for 
F  <  0.4. 

In  order  to  obtain  solutions  for  0.3  <  F  <  0.4,  we  used  scheme  III.  We 
found  that  scheme  III  gives  accurate  solutions  for  all  values  of  F.  As  n 
increases  the  coefficients  decrease  rapidly.  For  example 

d10  -  3.10"3,  d3Q  -  6.10*4,  d6Q  -  2.10-4,  dg0  -  2. 10~5  for  F  -  0.374.  In 
addition  the  solutions  of  scheme  III  agree  with  those  of  scheme  I  for  F  <  0.3 
and  with  those  of  scheme  II  for  F  >  0.4. 

In  Fig.  4  we  present  values  of  the  velocity  ~  at  the  separation 
point  S  versus  F.  The  velocity  is  equal  to  zero  for  F  <  0.36  and 

different  from  zero  for  F  >  0.36.  Therefore 

Fc  -  0.36  .  (20) 

This  value  is  about  40%  higher  than  the  value  guessed  by  Garabedian  . 

Profiles  of  the  jet  for  various  values  of  F  are  presented  in  Fig.  5. 

For  F  ■  •,  the  free  surface  reduces  to  the  vertical  line  y  ■  0.  As  F 
decreases  from  infinity,  the  jet  becomes  thinner.  As  F  approaches  zero,  the 
thickness  of  the  jet  tends  to  zero  and  the  free  surface  approaches  the 
horizontal  line  x  *  0.  For  F  >  F_  the  slope  of  the  free  surface  at 
x  ■  y  ■  0  is  vertical.  For  F  <  Fc  the  slope  at  x  -  y  -  0  is  horizontal. 

For  all  values  of  F,  the  free  surface  approaches  the  vertical  line  y  ■ 
j  *•  x  ♦  Therefore  the  jet  is  described  far  downstream  by  the  slender 

jet  theory  of  Keller  and  Geer3. 


This  theory  shows  that 


Values  of  the  dimensionless  velocity  —  at  the  separation  points  S  and 


u  -  iv  ~  -V(x)  as  x  ♦  -«* 


(21) 


Here  V(x)  is  real  and  positive.  Conservation  of  mass  and  equation  (4)  yield 
the  relations 

(1-2y)V(x)  -  1  (22) 

V2(x)  ~  -  2-  x  as  x  ♦  .  (23) 

F 

Eliminating  V(x)  between  (22)  and  (23)  gives 

x  - - - - r  .  (24) 

2<1-2y) 

Formula  (24)  is  the  asymptotic  shape  of  the  jet  far  downstream.  Our  numerical 

results  were  found  to  agree  with  (24)  for  x  large  and  negative.  This 

constitutes  an  important  check  on  our  numerical  scheme. 

The  solutions  of  Fig.  5  for  F  <  Fc  are  also  solutions  of  the  flow 

configuration  of  Fig.  1(a).  Therefore  there  exists  a  unique  "rising  bubble" 

for  each  value  of  F  smaller  than  Fc.  This  confirms  Garabedian ' s2  findings. 
2 

Garabedian  used  a  criterion  of  stability  to  suggest  that  the  unique 
physically  significant  "rising  bubble"  is  the  one  for  which  F  *  Fc.  Collins4 
reported  the  experimental  values  F  ■  0.25.  The  discrepancy  between  this 
experimental  data  and  our  theoretical  value  Fc  »  0.36  is  presumably  due  to 
the  three-dimensionality  of  the  real  flow  and  to  the  effect  of  surface 
tension.  The  bubble  profile  corresponding  to  F  ■  F  is  shown  in  Fig.  6. 


%>*  * 
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The  shape  of  a  two-dimensional  bubble  rising  at  a  constant  velocity  U  in  a 
tube  of  width  h  is  computed.  The  flow  is  assumed  to  be  inviscid  and  incom¬ 
pressible.  The  problem  is  solved  numerically  by  collocation.  The  results  confirm 
Garabedian ' s '2 )  findings.  There  exists  a  unique  solution  for  each  value  of  the 
Froude  number  F  =  U/(gh)l/2  smaller  than  a  critical  value  Fc.  Here  g  denote 


smaller  than  a  critical  value 


denote 


the  acceleration  of  gravity.  It  is  found  that  Fc  =  0.36.  In  addition  the  problez 
of  a  jet  emerging  from  a  vertical  nozzle  is  considered.  It  is  shown  that  the 
slope  of  the  free  surface  at  the  separation  points  is  horizontal  for  F  <  Fc 
and  vertical  for  F  >  Fc .  Graphs  and  tables  of  the  results  are  included. 
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